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The Asteroseismology of ZZ Ceti star GD1212 

Lin Guifang 1,2 , Li Yan 1,2 and Su Jie 1 2 3 

ABSTRACT 

The ZZ Ceti star GD 1212 was detected to have 19 independent modes from 
the two-wheel-controlled Kepler Spacecraft in 2014. By asymptotic analysis, we 
identify most of pulsation modes. We find out two set of complete triplets, and 
four sets of doublet which are interpreted as rotation modes with 1 = 1. For the 
other five modes, the four modes /i 3 , /i 5 , /i 6 and / 4 are identified as ones with 
l = 2; and the mode // is identified to be the one with 1 = 1. Meanwhile we derive 
a mean rotation period of 6.65 ± 0.21 h for GD 1212 according to the rotation 
splitting. Using the method of matching the observed periods to theoretical ones, 
we obtain the best-fitting model with the four parameters as M*/M 0 = 0.775, 
T e ff = 11400 K, log(Mn/M*) = —5.0, log(Mne/M*) = —2.5 for GD 1212. We 
find that due to the gradient of C/O abundance in the interior of white dwarf, 
some modes can not propagate to the stellar interior, which leads to the period 
spacing of the adjacent modes to become large. This feature is just proven by 
the observational data from GD 1212. All of these imply that GD 1212 may be 
evolved from an intermediate mass star. 

Subject headings: asteroseismology-stars:individual(GD1212)-white dwarf 


1. Introduction 


It’s well known tha t ab out 98 perce nt of stars will become white dwarfs in the end 
of their evolution (jWinget fe Kepler! 2008). Because there is no thermonuclear burning in 
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the interior of white dwarfs, the evolution of them is basically dominated by cooling. In the 
cooling curve, pulsating white dwarfs lie in three strips in the H-R diagram (jWinget & Kepler 
2008l h DOVs(75000K < T eff < 17000K), DBVs (22000K < T eff < 29000K), and DAVs 


DAVs or ZZ Ceti stars with the hydrogen atmosphere ( 

Bischoff-Kim & Metcalfe 

2011) 

to now, there are 148 ZZ Ceti stars ( 

Castanheira et al. 

2010). 


GD 1212 was confirmed to be a DA-type white dwarf (Holberg et al. 2002; Kawka 
et al. 2004), and was discovered to pulsate by Gianninas et ah (2006) with the main 
period of 1060.7s. In recent years, a series of observations by the spectroscopy and multi¬ 
color photometry were performed on GD 1212, and its atmospheric parameters(Kawka et 
al. 2004; Gianninas et al. 2006;Kawka et al. 2007; Subasavage et al. 2009;Sion et al. 2009; 
Gianninas et al. 2011; Giammichele et al. 2012) can be obtained. Hermes et al. (2014) did 
the photometric observations on GD 1212 using a two-wheel-controlled Kepler Spacecraft, 
and showed 19 independent pulsation modes ranging from 828.2s — 1220.8s. In their work, 
they only gave the discussions about the observed frequencies and various interpretation of 
the period-spacing, and did not match the pulsation modes with theoretical models. 


In this paper, we use the observed data of GD 1212 from the two-wheel-controlled 
Kepler Spacecraft (Hermes et al. 2014), and try to make the preliminary identification of 
19 independent modes in Section 2. In Section 3, we give the calculation of the theoretical 
model. By comparing the calculated periods of models with the observed ones, we give our 
best-fitting model. At last, the summary and discussions are given. 


2. preliminary identification of pulsation modes for GD 1212 


Mode identification is very important in the field of asteroseismology. In general, the 
pulsations observed in white dwarfs are non-radial g-mode ones. The pulsation mode can be 
characterized by three indices l,n, m, in which l is the spherical harmonic degree, n is the 
number of the radial nodes, and m is the azimuthal order. If a star is spherical symmetric, 
m would be degenerate, i.e. modes with the same l and n but different m will have the same 
frequency. If a star is rotating, rotation will destroy this spherical symmetry and lead to the 
frequency splitting. In the asymptotic approximation, the formula of frequency splitting due 
to rotation is as follows (Brickhill 1975) 


1Vl,n,m IVl,n + VTl\\ 


1 ( 1 + 1 ) 


]fi 


( 1 ) 


in which Wi^ m is the frequency of mode with indices l, n, m, is the rotation frequency of 
white dwarf. It is evident that a frequency with l = 1 will be splitted as three ones and called 
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a triplet; and a frequency with l = 2 will be splitted as five ones and called a quintuplet. 


According to the asymptotic theory of g-mode, on the other hand, an expression of 
periods of a mode with l and n(Unno et al. 1979; Tassoul 1980) is 

2vr 2 n 

Pl,n ~ 


\/l{l + 1) Jo r 


R N , 
— dr) 1 


( 2 ) 


in which N is the Brunt-Vaisala frequency, and R is the stellar radius. For a white dwarf, the 
right integral in the above equation is approximately invariable. This also means that pulsa¬ 
tion periods of different modes with the same l and two adjacent n should have an uniform pe¬ 
riod spacing. According to Eq.(2), there exists an approximate relation AP(1)/AP(2) rs./ V3 
between the period spaces of different modes with l = 1 and with 1 = 2. 


Using the above asymptotic expressions of g-mode, we reanalyze the 19 independent 
pulsation modes of GD 1212 from the two-wheel-controlled Kepler Spacecraft (see Table 1 
in Hermes et al. 2014 ). Owning to the frequency splitting of rotation, we identify 2 sets 
of complete triplets and 4 sets of doublet in Table [I] For 4 sets of doublet, we assume that 
the frequency with larger amplitude is the mode with m = 0. From Table [Q we can obtain 
that the frequency splitting of modes with l = 1 due to rotation is about 20 /iHz. According 
to Eq. (1) and each values of frequency splitting in Table [I] we can estimate the rotation 
period of GD 1212 which are listed in the last column in Table [0 In the calculated results, 
we can notice a rotation period of 8.08 h which is significantly different from others. This 
is because of the too small difference of the frequency between the modes of /j and /i 8 . 
According to Table 1 in Hermes et al. (2014), we notice that the error of frequency of /j 8 is 
3.54 /iHz. The small difference between f\ and /ig may be due to the observational error. So 
we assume that the mode of / 18 with frequency 857.14 /iHz and the other two modes of fu 
and fi make up of a set of triplet. At last, we estimate the mean rotation period of GD 1212 
is 6.65 ± 0.21 h, in which we don’t consider the rotation period by the frequency splitting 
between fi and /i 8 . 

In 19 pulsation modes, the other five frequencies are unidentified for the absence of 
frequency splitting. We perform asymptotic analysis for the five pulsation modes. From 
Table [2j we can get the period spacing of modes with l = 1 is about 37 s. According to 
the asymptotic theory of the non-radial g-mode , the period spacing of modes with l = 2 
is about 21 s. For the five unidentified modes, we notice that the period spacing of modes 
/i 3 ;/i 5 and fw is about 21 s, which can be identified as ones with 1 = 2. Because only one 
component for these modes has been observed, we hypothesize that they are the ones with 
m = 0. The period spacing of the last two modes /4 and /V is about 33 s, we cann’t firmly 
identify the two modes as ones with l = 1 or l = 2. In the following theoretical calculation, 
we fit them with / = 1 or l = 2. 
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Table 1: The possible frequency splitting for 1 = 1. ID is the frequency consistent with the 
one from Table 1 in Hermes et al. (2014). Freq is the frequency in unit of /tHz, d/ is the 
difference of adjacent frequency; Period is the observed period of modes in unit of s; m is 
the azimuth order, P rot is the rotation period of white dwarf in unit of hours derived by the 
frequency splitting. 


ID 

Freq(/iHz) 

d/(/tHz) 

Period(s) 

m 

Prot (h) 

/l4 

819.17 

20.79 

1220.75 

-1 

6.68 

h 

839.96 

17.18 

1190.53 

0 

8.08 

fis 

857.14 


1166.67 

+1 


fl7 

888.60 

21.85 

1125.37 

-1 

6.36 

h 

910.45 


1098.36 

0 


/l2 

920.71 

19.95 

1086.12 

-1 

6.96 

fb 

940.66 


1063.08 

0 


/20 

954.02 

21.29 

1048.19 

-1 

6.52 

f$ 

975.31 


1025.31 

0 


fn 

991.99 

21.18 

1008.07 

-1 

6.56 

h 

1013.17 


987.00 

0 


fs 

1166.17 

20.13 

857.51 

-1 

6.90 

fw 

1186.30 

21.15 

842.96 

0 

6.57 

h 

1207.45 


828.19 

+1 
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Table 2: The 11 observed frequencies with m = 0. ID is the frequency consistent with the 
one from Table 1 in Hermes et al. (2014). Freq is the frequency in unit of ^uHz, Period is 
the observed period of modes in unit of s; dP is the difference of observed adjacent periods; 
l is the spherical harmonic degree. 


ID 

Freq(/iHz) 

Period(s) 

dP(s) 

l 

h 

839.96 

1190.53 

92.17 

1 

/ 2 

910.45 

1098.36 

35.28 

1 

h 

940.66 

1063.08 

37.77 

1 

h 

975.31 

1025.31 

38.31 

1 

fz 

1013.17 

987.00 

144.04 

1 

fio 

1186.30 

842.96 


1 

h 

847.17 

1180.40 

33.28 

? 

h 

871.75 

1147.12 

190.25 

? 

fie 

1045.08 

956.87 

85.81 

2 

fl5 

1148.02 

871.06 

21.93 

2 

fl3 

1177.68 

849.13 


2 
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3. model fitting for GD1212 

For research of asteroseismology, it is essential to match the calculated periods of theo¬ 
retical models to the observed ones in order to find the best-fitting model. After the prelimi¬ 
nary identification of all the modes for GD 1212, we have a chance to do the astroseismology 
study for it. 


3.1. theoretical tools and input physics 

Firstly, we use the suite of open source MESA (Paxton et al. 2011, 2013) to get initial 
models of white dwarfs. Using the code of MESA, we adopt the metallicity Z = 0.02 and 
calculate the evolutions of stars with various masses from zero-age main sequence, to red 
giant branch, and asymptotic giant branch stages, until they evolve to become hot white 
dwarfs. The core composition prohles of the hot white dwarfs are the results of the stellar 
evolution before they become white dwarfs. 

Using the initial hot white dwarf models and incorporating the element diffusion (Thoul 
et al. 1994) into the White Dwarf Evolution Code(WDEC), we calculate the theoretical 
models of white dwarfs. WDEC was first introduced by Martin Schwarzschild (Schwarzschild 
& Harm 1965), and developed by other researchers(e.g. Kutter & Savedoff 1969; Lamb & 
van Horn 1975; Wood 1990; Montgomery 1998; Metcalfe 2001). The code adopts the tables 
of equation of state (EOS) from Lamb (1974) and Saumon et al. (1995) which are used in 
the degenerate core and partially ionized envelope, respectively. And the updated OPAL 
opacities of Iglesias & Rogers (1996) are used. In the theoretical calculations, we use the 
mixing length theory (Bohm & Cassinclli 1971) to treat the convection region, and adopt 
the mixing length parameter a as 0.6(Bergeron et al. 1995) . 

In order to calculate the models of DAV white dwarf, WDEC needs to input four 
parameters: the mass of white dwarf in unit of solar mass M*/M 0 , the effective temperature 
T e ff, the logarithmic hydrogen mass fraction log(Mn/M*), and the logarithmic helium mass 
fraction log(MH e /M*). The corresponding range of the four parameters are: 0.5 ^ M*/M 0 ^ 
0.90 with the step of AM*/M 0 = 0.005, 10850K ^ T c g- ^ 12250K with AT e g = 50K, 
— 10 ^ log(Mn/M*) ^ —4 with Alog(Mn/M*) = 0.5, and —4.0 ^ log(Mn e /M *) ^ —2.0 
with A log( Mn e /Iff*) = 0.5. According to the range of the four parameters, we notice that 
the four parameters cover the most part of the DAV white dwarf. With these grids, and 
taking the core composition prohles from MESA into account, we use WDEC to calculate 
all the models of white dwarf. Then we numerically solve the fully equations of linear and 
adiabatic oscillation in order to find each eigen-mode in theory by scanning in period . 
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3.2. the best-fitting model for GD 1212 


In order to find the best-fitting model for GD 1212, we compare the calculated periods 


for each model with 11 observed ones with previously identified spherical harmonic degree l. 
We introduce the expression of y 2 : 



( 3 ) 


where P c ^ is the calculated pulsation period from a model; P 0)i is the observed period; i is 
an integer which is from 1 to N, and N is the total number of observed periods. 

Using the above four-parameter space and the corresponding steps in Section 3.1, we 
calculate all the models of the white dwarfs and the pulsation periods from every model. 
Comparing the observed pulsation periods with the theoretical ones, the value of y 2 for 
each model can be calculated. We choose the model with the minimum y 2 as the best¬ 
fitting model, whose parameters are M*/M 0 = 0.775, T eS = 11400 K, log(M H /M*) = —5.0, 
log(Mn e /M*) = —2.5, and the corresponding minimum y 2 = 4.26. 

Figure 1 gives the grey-scale map of y 2 for the model with M* = O.775M 0 and \og(Mn e /Mfa 
—2.5. The abscissa is the effective temperature T eS , and the ordinate represents the loga¬ 
rithmic hydrogen mass fraction log(Mn/M*). From Fig. 1, we can see that the space of 
parameters with lower y 2 lies in the region with 10850K < T e ff < 11600K, and —6.0 < 
log(M H /M*) < -4.0. 

Because the calculated periods of the models can also be influenced by H/He mass 
fraction, we shorten the steps of the four parameters and reduce their range. The refined 
range of four parameters are as follows: 0.77 ^ M*/M 0 ^ 0.78 with AM*/M 0 = 0.005, 
10850K ^ T e ff ^ 11650K with AT e ff = 50K, —6 ^ log(Mn/M*) ^ —4 with A log(Mn/M*) = 
0.1, and —3.0 ^ log(M He /M*) ^ —2.0 with A log(M He /M*) = 0.1. We recalculate all the 
models and the corresponding pulsation periods. Finally we find that the best-fitting model 
is the same as the previous one obtained before. 

Table [3]gives the calculated pulsation periods of the best-fitting model. From this table, 
we notice that the unidentified mode fa with period 1147.12 s can be identified as a mode 
with l — 1, n — 27; the other mode j\ with period 1180.4 s can be identified as a mode 
with / = 2, n = 48. It can be noticed that the period difference between theoretical one and 
observational one for mode /A with period 956.87 s is 3.02 s and for mode f± with period 
1180.40 s is 4.08 s, which are comparatively large. However, the observational errors for 
these two modes are respectively 4.91 s and 4.02 s(Hermes et al. 2014). Our fitting errors 
are still comparable with the observational errors. 



11000 11200 11400 11600 11800 12000 12200 


Teff(K) 


Fig. 1.— The grey-scale map of x 2 f° r the model with M* = O.775M 0 and log(Mn e /M*) = 
—2.5. The abscissa is the effective temperature T e g, the ordinate represents the logarithmic 
hydrogen mass fraction log(M H /M*). 
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Table 3: The periods of the best-fitting model. Data of the first and fifth columns are the 
values of / and n for the calculated periods. P c is the calculated pulsation period; P Q is the 
observed period. The difference | P c — P Q | of them is also given. 


/ = 1 

n 

PcW 

P.(s) 

(ia-Poi) 

/ = 2 

n 

Pc( 8) 

P.W 

(1 Pc -Pol) 

15 

694.12 



34 

827.75 



16 

722.13 



35 

849.31 

849.13 

0.18 

17 

761.38 



36 

871.62 

871.06 

0.56 

18 

795.63 



37 

897.73 



19 

842.97 

842.96 

0.01 

38 

924.59 



20 

884.24 



39 

953.85 

956.87 

3.02 

21 

919.95 



40 

982.76 



22 

943.47 



41 

1005.79 



23 

986.57 

987.00 

0.43 

42 

1016.38 



24 

1022.71 

1025.31 

2.60 

43 

1042.68 



25 

1060.40 

1063.08 

2.68 

44 

1069.46 



26 

1099.05 

1098.36 

0.69 

45 

1096.79 



27 

1148.18 

1147.12 

1.06 

46 

1123.37 



28 

1192.78 

1190.53 

2.25 

47 

1150.43 



29 

1232.93 



48 

1176.32 

1180.40 

4.08 

30 

1256.67 



49 

1200.64 
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3.3. Constraints from the theoretical models 

Figure 2 gives the period spacing of calculated consecutive modes with / = 1 and l = 2. 
The upper panel is the one with / = 1 and the lower panel is the one with l = 2. The hollow 
circles represent the period spacing of theoretical modes; and the filled circles denote the 
one of observed periods with adjacent modes. For the panel with l = 1, the number of radial 
nodes n is from 19 to 31, and for the other with / = 2, n is from 33 to 53. In this figure, 
we notice that results of theoretical models are in good agreement with the observations. 
Furthermore, we notice that the first three observed period spacings are around 37 s, then 
the next two observed period spacings suddenly increase to 48.76 s and 43.41 s. Such large 
increase in the period spacing usually implies changes of propagation region for different 
pulsation modes. We will discuss this feature in the following. 

Figure 3 shows the composition profiles in the best-fitting model with M* = O.775M 0 . 
We can see that there are 2 transition zones in the interior of the white dwarf. One is the 
transition zone of H and He which lies in the position of log(l — M r /M*) ~ —5. The other is 
the one of He and C/O which is located at the position of log(l — M r /M*) ~ —2.5. It should 
be noticed that the core of the white dwarf is the result of nuclear burning and convective 
mixing in the stellar evolution. Before the star becomes a white dwarf, a convection core 
lies in the central part of the star, which makes the elements completely mix in the core. 
Outside the convection core, there is a region where a gradient of C/O abundance forms. 

In order to explain the feature that the period spacing of consecutive modes with / = 1 
and n = 26 suddenly becomes larger in Fig. 2, we calculate the so-called modal mass M mode , 
whose expression is: 

f R 

Mmode = 4vr / [| | 2 +/(/ + 1) I \ 2 }por 2 dr (4) 

Jo 

in which is the radial displacement, £h is the horizontal displacement, po is the stellar 
density, and r is the radius. 

The distributions of the so-call modal mass for some modes with / = 1 are shown in Fig. 
4. The abscissa is the nondimensional quantity log(l — M r /M*). The ordinate represents a 
normalized integrant of the so-call modal mass. The panels from top to bottom are the ones 
for modes with n from 23 to 30. The scales in the ordinate and abscissa are the same for 
all panels. From this figure, we can notice that most of pulsation energy for all modes are 
distributed in two transition zones. One of them lies in the position of log(l — M r /M*) ~ —5, 
which is the transition zone of H and He. The other is around log(l — M r /M*) ~ —2.5, which 
is the one of He and C/O. We also find that the pulsation energy in the interior of C/O core 
for modes with n = 26 and n = 27 is very small compared with others. Meanwhile their 



Period(s) 


Fig. 2.— The period spacing of consecutive modes with l = 1 (upper) and l = 2 (lower). 
The hollow circles represent the period spacing of theoretical modes and the filled circles 
denote the one of observed periods with adjacent modes. 
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period spacings with respect to the adjacent mode increase suddenly. This feature is evident 
very clearly also in the observed data (see Fig. 2). On the other hand, the pulsation energy 
in the core for the mode with n = 29 is very large, and the corresponding period spacing 
is very small. All of these can be explained by Eq. (2). When the pulsation energy can be 
propagated into the stellar interior, the range of the integral in Eq. (2) must be extended to 
the stellar core. The integral part becomes very large, which leads to a small period spacing, 
and vice versa. 

Comparing Fig. 3 with Fig. 4, we notice that the region with such a gradient of C/O 
abundance lies in the position of log(l — M r /M*) ~ —0.5. It separates the stellar core into 
two parts: the inner core has a lower C/O ratio and the outer core has a higher C/O ratio. 
Due to the gradient of C/O abundance, some modes can not propagate to the stellar interior. 
This leads to the period spacing for these modes to become large. This feature is just proven 
by the observational data from GD 1212. It means that the gradient of C/O abundance 
should be exist in the interior of GD 1212. For GD 1212, the mass of the best-fitting model 
is M*/M 0 = 0.775. In the central He burning phase, there is an convection core where the 
elements are to be mixed completely. The profiles of C/O elements outside the convection 
core will be formed by the shell burning of He element in the phase of AGB. All of these 
imply that GD 1212 may be evolved from an intermediate mass star. 


4. Summary and discuss 

Asteroseismology provides a unique window to detect the internal structure of the pul¬ 
sating stars by matching the observed periods to theoretical periods generated by adiabatic 
pulsation. 

In this paper, by asymptotic analysis, we identify most of pulsation modes of GD 1212 
observed by the two-wheel-controlled Kepler Spacecraft (Hermes et al. 2014). Because of 
the rotation splitting, we identify two set of complete triplets and four sets of doublet. We 
interpret them as rotation modes with / = 1 owing to the resemble frequency splitting. Thus 
we estimate a mean rotation period is 6.65 ± 0.21 h for GD 1212. Using asymptotic analysis 
of non-radial g-mode, we identify the three modes / 13 , /15 and /i 6 as ones with 1 = 2. The 
two modes / 4 and f 7 are identified as ones with / = 2 and 1 = 1, respectively, by matching 
the observed periods with theoretical ones. By the asteroseismology study, four parameters 
of the best-fitting model for GD 1212 can be obtained as M*/M 0 = 0.775, T e s = 11400 K, 
log(Mn/M*) = —5.0, log(Mn e /M*) = —2.5. Simultaneously, we find that due to the gradient 
of chemical element abundance outside the convection zone, some modes can not propagate 
to the stellar interior. This leads to the period spacing of the adjacent modes to become 
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Fig. 4.— The distribution of the fraction of the so-call modal mass. The abscissa is the 
nondimensional quantity log(l — M r /M*). The ordinate represents a fraction of the so-call 
modal mass. The panels from top to bottom are the ones for modes with / = 1 and n from 
23 to 30. 
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large, which is just proven by the observational data from GD 1212. This means that the 
gradient of C/O abundance should exist. And ZZ Ceti star GD 1212 may be evolved from 
an intermediate mass star. 

In order to compare the theoretical parameters with ones obtained from observation for 
GD 1212, we collect the values of the parameters from the various observations. In Table 
|4j we list the values of parameters: the effective temperature T eff , logg, and the mass M* 
obtained from the spectrum and multi-color photometry. The last row is our calculated 
results for GD 1212. From this table, our result is very close to the ones from Subasavage 
et al. (2009) and Giammichelc et al. (2012).In future, further observations of GD 1212 are 
needed. This will help us to find more independent modes and then further constrain the 
model of GD 1212. 
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Table 4. 


The values of the parameters from the various observations. 


Teff(K) 

log g 

M/M© 

log M h /M* 

log M He /M* 

Ref 

10960T75 

8 .2T0.1 

0.73T0.06 



1 

11040T132 

8.11T0.038 

0.67 ±0.02 



2 

11010 T210 

8.05T0.15 

0.63T0.09 



3 

11000T300 

8.25T0.04 

0.76T0.02 



4 

11270T165 

8.18T0.05 

0.71T0.03 



5 

10938T317 

8.25T0.03 

0.76T0.02 



6 

11400 

8.37 

0.775 

-5.0 

-2.5 

7 
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